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We use a version of the density functional theory to study the changes in the height of the tethered layer 
of chains built of jointed spherical segments with the change of the length and surface density of chains. 
For the model in which the interactions between segments and solvent molecules are the same as between 
solvent molecules we have discovered two effects that have not been observed in previous studies. Under 
certain conditions and for low surface concentrations of the chains, the height of the pinned layer may attain 
a minimum. Moreover, for some systems we observe that when the temperature increases, the height of the 
layer of chains may decrease. 
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1. Introduction 

The materials modified with tethered chains play an important role in several technological 
fields. With the proper choice of the grafting surface, the length of chains and the grafting den- 
sity, one can obtain novel products exhibiting desired behavior. The materials modified in such 
a way are utilized in a wide range of applications including chromatography, adhesion, lubrica- 
tion and colloidal stabilization |TH3] ■ Recent studies have focused on the development of "smart" 
materials, whose properties can be readily altered by applied stimuli. They can find applications 
in drug delivery, as chemical sensors and in controlling the nanoparticle transport [4-7 . Further- 
more, applications in nanotechnology include end grafted polymers on nanopatterned surfaces, on 
nanoparticles, or on carbon nanotubes [SHU]- 

Much effort has been directed toward theoretical description of grafted polymer layers. The 
studies included scaling theories [IS'-TS^, classical self-consistent field methods [17, T? , single chain 
mean-field methods LlSj, density functional theories |18ff2?[ and computer simulations [23:-34:. . In 
our opinion, particularly important for further theoretical studies of systems involving tethered 
chains is the development of density functional theory (DFT) fIU\ [2TJ [55H1^ that is based on the 
approach proposed by Yu and Wu [43 . In general, this type of theory represents a powerful tool for 
the study of adsorption and phase behavior of confined systems. Within this approach, theoretical 
description of "pinning" the terminating segments of chain molecules to the walls can be realized 
by using the ideas of the previously developed methods of bonding the fluid species to the walls in 
an associative treatment |44l |45] . 

Recent theoretical works based on the DFT |20| have mainly concentrated on the 

description of adsorption of fluids on modifled surfaces, while the problem of the changes in the 
layer of pinned chains itself has attracted less attention. Of course, the accumulation of fluid 
molecules near the adsorbing surface induces changes in the tethered layer, especially when the 

* Dedicated to Professor Yu. Kalyuzhnyi on the occasion of his 60th birthday. 
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adsorbed fluid undergoes prewetting of layering transitions [101, However, no systematic density 
functional study of changes in the structure of tethered chains has been undertaken so far. 

The problem of a change of the conflguration of tethered chains with their length, surface density 
and the quality of a solvent has been considered in a great number of publications. One of the first 
theoretical works devoted to that problem were the publications by Alexander and de Gennes |131 
I14| . However, the references to numerous more recent works can be found in the review article of 
Descas et al. [46| and in the PhD thesis by Chremos |47|. We know that in a good solvent, for not 
too high surface densities of tethered chains and when the free segments (i.e. all but the tethered 
segment) are not too strongly attracted by the surface, the chains assume the so-called mushroom 
configuration. With an increase of the surface density, the overlapping mushroom configuration is 
observed. This configuration precedes the brush regime that occurs at high surface concentrations. 
Despite the great number of experimental, theoretical and computational investigations on the 
behavior of end anchored chains, the quantitative determination of the parameter space appropriate 
for detailed studies of transitions between possible configurations of tethered chains is still vague. In 
addition, the predictions of theories do not always coincide with the results of simulations |32lH5] . 
Moreover, in our opinion insufficient attention has been paid to some problems, such as the effect 
of competition between strong adsorption of fiuid molecules and adsorption of free segments of a 
chain on the structure of the chains and to the effect of the architecture of branched polymers on 
the structure of the formed layer |49) . 

The description of systems comprising tethered chains at a microscopic level requires the knowl- 
edge of numerous parameters. One has to know not only the surface density and the length of the 
chains, "quality" (and density) of the solvent and the temperature, but also details about the archi- 
tecture of the chains, parameters describing the interactions between all the components themselves 
and between the components and the surface, and, in some cases, the topography of the end an- 
chored chains on the surface. A huge number of parameters of the model delimits the possibility of 
an effective scanning of the space of the model parameters for a precise evaluation of the regions 
of different configurations of chains and the transitions between them by long lasting computer 
simulations. Such studies, however, can be carried out by using theoretical approaches based on 
microscopic models. Therefore, the aim of the present work is to employ the DFT used in our 
previously works |20| |2I] to investigate how the thickness of the pinned layer changes with 

their length and with their surface density. Our principal purpose is to check whether the theory 
employed by us is capable of reproducing the results of coarse-grained approaches |13fll5| . We are 
also looking for some new effects that have not been observed in previous studies. 

2. Theory 

Following our previous works |20l I^TJ [321 HH we consider a fluid of spherical molecules of 
diameter cTf in contact with a modified substrate. The surface is covered by the film of preadsorbed 
chain molecules. All the chains are represented by tangentially jointed M spherical beads of the 
same diameter, a^- The chain connectivity is enforced by the bonding potential between nearest- 
neighbor segments, Vh- This potential satisfies the equation [13] 

exp[-/3V^b(R)] = n ^d'^^ ~ - ^c)/4^'Te (1) 

2 = 1 

where R = (ri, r2, . . . , tm) is a vector specifying segment positions of consecutive segments. Each 
polymer molecule of the chemically bonded phase contains one surface-binding segment located at 
its end (indexed as "1") that interacts with the wall via the potential 

eM-l3vsi{z)]=Cd{z-aj2) (2) 

where z is a distance from the surface and C is a constant. This potential implies that the surface- 
binding segments always lie at the distance x = (7c/2 from the surface. The remaining segments of 
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the grafted molecules (j = 2, 3, . . . , M) are "neutral" with respect to the surface, i.e., they interact 
with the surface via a hard-wall potential 

The fluid molecules, however, interact with the surface via the Lennard-Jones (9,3) potential of 
the form. 

v{z) = £fe [{z/zof - {z/zof] , (4) 

with zq — (7f/2. 

We assume the Lennard-Jones (12,6) type interactions between all segments, between the seg- 
ments and fluid particles and between the fluid particles 

u,j{r) = e,j[{r/a.,jr-{r/a,,)'], (5) 

where (Tij = 0.5((Ti -I- crj), i,j = c,f. Employing the perturbative treatment we split these inter- 
actions into the repulsive (reference) and attractive (perturbation) parts according to the Weeks- 
Chandler- Anderson scheme |50| Uij{r) — Uij,-c[{r) + u.y.pcrt('")i where 

In order to proceed, let us introduce the notation, /9*^'^)(R) and p{z), for the density distribution 
of chains and of spherical species, respectively. However, the theory is constructed in terms of the 
density of particular segments of chains, Psj{z), and the total segment density of chains, Ps{z). These 
densities are introduced via commonly used relations, see the original development in |20 [ l2l ]l39l - l42| 

M M „ 

p.(r) - 5]ps,(r) = dR<5(r, - r)p(^)(R). (7) 

In the system under study all the local densities are the functions of the distance from the surface 
only. 

The system is studied in a grand canonical ensemble with the constraint of constancy of the 
number of chain molecules, i.e. 

Pc = J dzpsj{z) (8) 

where pc is the number of chain molecules per area of the surface, i.e. it is the surface density of 
the chains. Of course, the integral in equation Q does not depend on the segment index j. 
The thermodynamic potential appropriate to the description of the system is 

Y = F[p,{z),piz)] + / dzpizMz) -p]+J2 

dzpsjVsj{z), (9) 

where F[ps{z), p(z)] is the Helmholtz free energy functional and p is the chemical potential of the 
reference system fluid. The expression for F[ps{z) , p{z)] is obtained from the theory described in 
our earlier works |20 [ [2T | 139^42] and in the original papers by Yu and Wu [43]. The density proflle of 
fluid molecules and the segment density proflles are obtained by minimizing the functional (|9| under 
the constraint ([s]). For the sake of brevity we do not present the resulting density profile equations. 
They can be found in the above cited publications. We should only note that the reference system 
comprises only fluid species. The chemical potential of the fluid, p is thus given by 

p^ kT\npk + Phs + Pb J "fF,pert(?')dr, (10) 
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where ph is the reference system density and /ihg is the chemical potential of hard-spheres of 
diameter Cf , resulting from the Carnahan-Starling equation of state. 

The system under study is described by several parameters. In order to reduce the number 
of parameters to a minimum, we assume that — ac = (Td = cr and that — £cf = £■ 
The Lennard-Jones (12,6) potential parameter, characterizing the interactions between segments 
of the chain, £cci as well as the fluid-solid energy parameter are expressed in the units of e, 
£*j. = £cc/e and £fj, = £fs/£. The reduced densities are defined as usual, p'^ = P\^cf^, P*{z) = p{^)'^^i 
plj{z) = psj{z)a'^ , P*s{z) = Ps{z)(j'^ and p* = pccr'^ ■ Moreover, we introduce the reduced temperature 
as T* = kT/e. 

Almost all the calculations (unless otherwise explicitly stated) have been carried out for the fluid 
bulk density = 0.66. This density nearly coincides with the liquid density on the liquid-vapor 
coexistence curve at T* = 1. 



3. Results and discussion 

Before discussing the results we should emphasize that it is unclear how to determine the brush 
height in a precise meaning from the density profiles obtained via DFT or from simulations. It 
has been a widespread practice to define the brush height, {h), from the first moment of the total 
segment density profile |26} [5l] [52] 

^dzzp^jz) 

{h)^a-jr- — . (11) 

Usually, (cf. [261 15T]) the coefficient a was set to be a = 8/3, but also a = 2 was used (see, e.g. |52j ) 
and the latter value has been adopted in our calculations. Moreover, usually the height of the 
brush, as well as the surface density of tethered chains were scaled by the radius of gyration (or by 
some powers of the radius of gyration) of the chain molecules, cf. |26l 15^ 1151 [511 [5^ . Unfortunately, 
the value of the radius of gyration cannot be self-consistently obtained within the framework of 
the considered approach. Therefore, in our calculations we use the reduced units as defined at the 
end of the previous section. 



3.1. Change of the brush height with the length of chains 

We begin with the discussion of the changes of the brush height with the length of the tethered 
chains. In figure [l] we show the plots of {h)/a versus M for a number of systems characterized by 
different values of the intermolecular potentials and obtained for different surface grafting densities. 
In part (a) the solid lines were obtained assuming weak interactions between ffuid molecules and 
solid surface, e^^ = Sis/e — 1. The remaining lines were evaluated for solid surfaces attracting 
the ffuid molecules more strongly, see the ffgure caption. The surface density of grafted chains is 
p* =0.1 for all, but the one curve, for which p* — 0.02 (the solid line decorated with solid circles). 
From first glance, for M > 10, almost all the curves plotted here are well approximated by straight 
lines (an exception is the solid line decorated with solid circles, calculated for p* = 0.02). However, 
much better approximations are obtained having assumed a power-law dependence, cf. |53| 

(h) cx M". (12) 

Even for the segment-segment interactions of a hard-sphere type (solid line decorated with stars) 
the approximation of numerical results by a straight line of the form (h) = a + bM yields the 
correlation coefficient Rc of the order of i?c ~ 0.990, while the use of equation ( [T2| ) gives Rc ~ 0.9999 
for a = 0.9134. 

When the attractive interactions between the segments are turned on, the exponent a decreases. 
For £cc = 0.6£ we obtain a = 0.850 (solid line with crosses; the correlation coefficient in this case 
is again of the order of 0.9999), while for Scc — e we have a = 0.738 (solid line with open squares; 
the correlation coefficient is now slightly lower, Rc ~ 0.999). An increase of the ffuid-solid energy 
decreases slightly the value of a. For example, when £cc = 0.6£ and e^^ = 12 (dashed line with 
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(a) (b) 

Figure 1. The dependence of {h)/a on M. In part (a) the nomenclature of the Hnes decorated 
with different symbols is as follows. Solid lines have been evaluated for = 1 and for the chains 
built of hard-spheres, e*^ = (stars); e*^^ — 0.6 (crosses) e*c = 1 (open squares). Dotted line 
with filled squares is for — 5 and elc ~ 1- Dashed lines are for ej^ — 12 and for e*c = 0.6 
(crosses) and e'^ ~ 1 (open squares), while dash-dotted line with triangles is for ej, = 18 and 
e*c = 1. The values of the remaining parameters are: T* = 1 and p* = 0.1. Additionally, solid 
line with filled circles has been obtained for — 1 and e^c = 1, T* = 1 and = 0.02. 
In part (b) we show the effect of surface density p* and of the temperature on the evaluated 
curves {h)/a = f{M). Solid lines are at T* = 1 and p* = 0.1 (open circles), p* — 0.12 (open 
squares), p* — 0.15 (open diamonds) and p* = 0.20 (open triangles). The dashed line with 
stars is for p* — 0.1 and at T* = 4. The values of the remaining parameters are Sf^ = 12 and 

£* = 1 



crosses) we have a — 0.834, instead oi a = 0.850 for ~ 1 (solid line with crosses). Similarly, for 
£cc = £ and for — 18 (dash-dotted line with triangles) we obtain a — 0.720. The last value is 
lower than that obtained for Scc = £ and Ef^ = 1. The effect of the parameter e^^ on the exponent 
a is much weaker that of the parameter Scc ■ 

The solid line decorated with filled circles that has been evaluated for much lower surface density 
of tethered chains deviates very much from all the remaining results. In this case the exponent of 
equation ( 12 1 is much lower and equals a w 0.49. This indicates that the surface density of tethered 
chains greatly effects the scaling relationship (12 1. To investigate this problem, we calculated the 
values of {h)/a for different values of p*, keeping the values of the following parameters constant: 
£cc — £ and e*^^ — 12, see figure [T| (b). Solid lines displayed here correspond to T* = 1, while the 
dashed line is for T* = 4. For p* = 0.1 (solid line with open circles) we found that a = 0.729, 
while for p* = 0.2 (solid linear with triangles) - a — 0.828. Therefore, we conclude that an increase 
of the surface density of tethered chains increases the value of a. This is a rather obvious result 
because in the case of higher values of p* the chains are more stretched and there is "less room" 
for lateral motions of the segments which could possibly lead to "coiling" the chains. However, the 
effect of the temperature on the exponent a is not so obvious. The dashed line was obtained at 
T* = 4 for p* =0.1. Now, a « 0.90. For lower values of M the height of the brush at T* = 4 is 
lower than at T* = 1, but for longer chains the situation is reverse. One can argue here that the 
average brush height, (h), should be related to the end-to-end distance, Ri-n, of a single polymer 
chain, at least at low grafting density. Unfortunately, the density functional approach used in this 
work does not allow for the evaluation of the end-to-end distance, and therefore we cannot carry 
out any test for a relation between (h) and Ri-n using self-consistent data in the framework of 
the applied theory. 



33604-5 



M. Borowko et. al 



The configuration assumed by the chains depends on entropic as well as on energetic contribu- 
tions to the free energy. The chains maximize their configurational entropy by adopting a random- 
walk-like configurations. However, depending on the energies of segment-segment and segment-fiuid 
attractions, they tend to maximize the contact between relevant species. Moreover, the attraction 
between fiuid molecules and the solid leads to accumulation of the fiuid molecules in the vicinity of 
the solid surface. At T* = 4 the Boltzmann factor of the fiuid-solid interactions is much lower than 
at T* = 1. Therefore, the accumulation of fluid molecules and thus the effect of "superseding" the 
segments from the region close to the wall is smaller at T* = 4. The adsorption of fluid molecules 
on the brush structure plays a more significant role for shorter chains. For longer chains, however, 
the entropic effects prevail and at T* = 4 the brush becomes more stretched than at T* = 1. We 
return to the problem of temperature dependence of the brush height hereinbelow. 




(a) (b) 

Figure 2. Part (a). Total segment density profiles (upper panel), the density profiles of the last 
segment (middle panel) and the fiuid density profiles for the chains of different length. Dotted 
line is for Ad = 8, dash-dotted line - for AI = 16, dashed line is for AI — 24 and solid line is 
for AI — 32. Part (b) shows the density profiles of individual segments, plj{z) for the values of 
j given in the figure. Here the chain length is M = 32. The inset to part (b) shows the scaled 
segment density profiles, as described in the text. The consecutive curves from left to right are 
for j — 2, 6, 10, 14, 18, 22, 26 and 30. The segment-segment interactions are of hard-sphere 
type, p* = 0.15, T* — 1 and = 1. 

We have monitored the changes in the structures of the brush and of the fluid with the changes 
of the parameters of the model by inspecting the density profiles. Examples of the obtained results 
are shown in figures |2] and |3] Figure |2] (a) shows how the total density profiles (the upper panel), 
the profiles of the last segment (the middle panel) and the fluid density profiles (the lower panel) 
change with the length of tethered chains. We see that if the length of the chains increases, the 
initial parts (i.e. the part close to the wall) of the total segment density profiles and of the fluid 
density profiles remain almost the same as for shortened chains. An increase in the number of 
segments M causes expansion of the region where the segments and the fluid molecules exhibit a 
layered structure. Instantaneously, the "tails" of the profiles become more diffuse for longer chains. 
In the middle part of figure |2] (a) we also see that even for long chains there is nonzero probability 
of finding the last segments in the vicinity of the wall. Such a shape of the segment density profiles 
suggest coiling of the chains which can even lead to the formation of loops with both terminating 
segments of a chain located at the wall. 
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Figure |2] (b) displays the profiles of individual segments with even numbers for tethered 32- 
mers. With an increase of the segment number, the profiles become more and more diffuse. Ac- 
cording to classical self-consistent field theory, the re-scaled distributions of the segments p'j{z') = 
plj{z) sin[7rj72M] should be an universal function of z' = z/sin[7rj/2Af], i.e., the plot of p'^ versus 
z' should be independent of the segment index j [54H56] . The inset to figure [2] (b) provides a test 
whether the above scaling can be also applied to the DFT results. Evidently, the self-consistent 
field scaling fails when applied to the density functional results of this work. 




0.006 - 




Figure 3. Part (a). Total segment density profiles (upper panel) and fluid density profiles (lower 
panel) of M = 80-mers for = 5 (solid lines) and Efg = 15 (dotted lines). In part (b) we display 
the profiles of the last segment (upper panel) and of the middle segment (j — 40) and of the 
second segment, j = 2 (the latter are divided by 10). The values of the remaining parameters 
are RT = 0.1, T* = 1 and Ecc = £■ 



Figure [3] (a) illustrates the change of the total segment density profile of tethered chains and of 
the profile of the fluid (lower and upper panels, respectively), while figure|3](b) shows the changes of 
the last (lower panel) and of the middle (upper panel) with an increase of the fluid-solid attraction. 
We see that when s^^ increases, the total segment density proflle is pushed away from the solid 
surface by the adsorbing fluid. This effect is more pronounced for the last segment than for the 
initial one of middle segments of the chain (figure [S] (b)). 

3.2. The role of the surface density of the chains 

The next problem we want to study in detail is how the surface density of the pinned chains 
effects the height of tethered layer. We begin our discussion with the results for rather short chains 
built of 16-mers (the solids modified with the chains of a comparable length are widely used as col- 
umn fillings in chromatography). In figure |4] we demonstrate the effect of some selected parameters 
of the model on the relationship of (h) / a upon p* . Three curves displayed in part (a) have been 
calculated assuming that the segment-segment interactions are the same as the segment fluid and 
fluid-fluid interactions, £cc — £■ The temperature was kept constant and equal to T* = 1. The curve 
decorated with open squares is for the solid that weakly attracts fluid molecules, e^^ — l, while the 
curve decorated with open diamonds is for strongly attracting solid surface, = 15. These two 
curves were evaluated assuming that the fluid density was pj^ = 0.66. However, we also present here 
the results for tethered chains in a vacuum (curve decorated with solid circles) . Part (b) illustrates 
the effects of weakening of the segment segment interactions and of the temperature on {h)/a. 

For very low surface densities of tethered chains, the brush height remains almost constant. 
This is not surprising, because under such conditions a tethered chain does not almost "feel" the 



33604-7 



M. Borowko et. al 




(a) (b) 

Figure 4. The dependence of (h) /a on the surface density of tethered 16-mers. The nomenclature 
of the lines in part (a) is as follows: — 15 and = 0.66 - line with open diamonds; ej^ — 1 
and = 0-66 - line with open squares. The line with solid circles was obtained for the brush 
in a vacuum, p'^ = 0. In all cases Ecc = £ and T* — 1. Dotted straight line has the slope of 1/3 
and fits all the results in the limit of high Pc, p* G [0.5,0.9]. In part (b) the consecutive lines 
were obtained for the following parameters: Ecc ~ 0.6e and T* = 1 - line with stars; Ecc = 
(hard-sphere segment-segment interactions) and T* = 1 - line decorated with crosses; Ecc = £ 
and T* — 5 - line decorated with triangles. In all cases — 0.66 and Ef^ — 1. The dotted and 
dashed straight lines denote the approximations of the log-log plots, as described in the text. 



presence of other chains and assumes the configuration that is only effected by the presence of 
solvent molecules. The accumulation of solvent molecules in the vicinity of the solid surface leads 
to an increase of (h) /a. Of course, the last statement is true if the solvent-fluid attraction is similar 
(in our case identical) to the fluid-fluid attraction. The behavior of the brush in "bad" solvent can 
be different [M], but this problem requires additional studies. 

With further increase of p* and when the attractive fluid-solid potential is weak, there develops 
a minimum on the curve of {h)/a vs. p* (see the line with open squares in flgure |4] (a)). In this 
case, the minimum is located at p* « 0.1. Similarly, a minimum on the curve of {h)/a vs. p* is 
observed if the brush is in contact with a vacuum, but now it is located at much lower surface brush 
density, p* ~ 0.04. For a strongly attractive wall, this minimum disappears (solid line with crosses 
in figure |4] (a)). Also, the minimum vanishes when the segment-segment attraction is lowered or 
when the temperature is raised, see figure |4] (b). Evidently, this phenomenon is caused by the 
segment-segment attractive interactions. At low temperatures, the attraction between segments 
causes coiling of the chains near the wall. However, in the case of the wall strongly attracting fluid 
molecules, the coiling is inhibited by the fluid molecules accumulated at the wall. 

In general, two characteristic regions on the curves {h)/(T vs. p* in flgures|4](a) and[4](b) can 
be distinguished: the region of almost constant values of {h)/a at low surface densities of tethered 
chains and the region where the logarithm of the height of the brush changes with In p* almost 
linearly. Such a behavior was also observed in computer simulations |26l I48j . 

Classical scaling theories predict that the height of the brush scales with the surface brush 
density as 

{h) « {p,r (13) 

with 7 = 1/3, as predicted by the self-consistent mean-filed theory (see, e.g., (Ml), or with 7 w 0.35, 
as predicted by the Alexander-de Gennes approach (see, e.g., |48]. The dotted line plotted in 
figure [4] (a) has the slope of 1/3 and it well fits all the data for 0.5 < p* < 0.9. However, for the 
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values of p* higher than 0.9 the changes of {h) /a are much smaller than predicted by equation ( 13 I. 
Of course, for a very dense brush, its height attains a maximum corresponding to fully stretched 
chains. 

At high temperature, the linearity range of the log-log plot is much wider, cf. dotted line in 
figure |4] (b) . This line well approximates the {h{p*)/a) curve evaluated at T* = 5 (the line with 
triangles), but now the exponent 7 w 0.43 is significantly higher than predicted by the theory of 
Alexander and de Gennes |13fll5| . Similarly, when the segment-segment interactions are of hard- 
sphere type, there exist a range of surface densities p*, within which the log-log plot is linear 
(dashed line in figure |4] (b)), but now the slope is much lower, 7 w 0.14, than predicted by the 
classical theories. Only when Ecc — 0.6e (the curve with stars in figure |4] (b)), the straight line 
approximating the data in a relatively wide range of p* values has the slope close to 1/3 (the 
relevant line has been omitted in figure |4](b)). 




Figure 5. Examples of the total segment density profiles of 16-mers and the profiles of fluid 
calculated along the curves displayed in figure 4 (a). Part (a) is for the brush in contact with 
a vacuum (line with solid circles in figure 4). Upper panel shows the results for p* — 2 ■ 10~^ - 
dotted line, 0.04 (this surface density corresponds to the minimum of {h)/a) - dashed line, 
0.4 - dash-dotted linear and 0.7 - solid line. Lower panel is for = 0.9 - dashed line and for 
1.15 - solid line. Part (b) was calculated along the line with open squares in figure 4 (a). Upper 
panel displays the total segment density profiles and lower panel displays the profiles of fiuid. 
The calculations are for p* — 2 ■ - dotted lines, 0.1 (this surface density corresponds to 

the minimum of {h}/a) - dashed line, 0.4 - solid line and 0.9 - dash-dotted line. Part (c) was 
calculated along the line with open squares in figure 4 (a). The abbreviation of the lines and the 
surface densities of tethered chains are the same as in part (b). 

Figure[5]illustrates how the structure of the brush built of 16-mers in a vacuum changes with p*. 
At extremely low surface density, p* = 10^^, the brush extends to the distance around z/a w 5.5 
and the profiles pl{z)/ p* remain almost unchanged up to p* « 5 • 10~^ (the relevant results have 
been omitted for the sake of brevity) . Then, the coiling process starts due to attraction between the 
segments and the minimum height of the brush is attained at p* k, 0.04. The profile corresponding 
to the brush of minimal height is shown by the dashed line in the upper panel of figure [5] (a) . 
This curve demonstrates that the formed structure essentially comprises two layers at the wall. 
Further increase of p* leads to the brush expansion. Initially, the layer at the wall resembles a 
liquid-like film, but then a layered structure develops. For p* < 0.9, the number of layers is lower 
than the number of the segments M, but at very high values of p*, the number of layers increases 
up to 2M (see the lower panel in figure |5] (a)). The brush layer tends to assume the solid-like 
configuration |57j . Unfortunately, one-dimensional density functional theory used in this work fails 
in that region of the surface brush densities. 

The presence of the fluid changes the structure of the brush, cf. figures [5] (b) and|5](c). When 
p* is not too high (say, p* < 0.4, the fiuid itself forms layered structure, especially for higher values 
of Efg . The presence of a well-pronounced fiuid layers in the vicinity of the wall inhibits coiling of 
the brush and breaks the development of a minimum on {h{p*)/a) curve. The layers of the fluid 
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also cause an increase of the brush height, in comparison with the brush in a vacuum. However, for 
higher surface densities of tethered chains, the fluid molecules are no longer capable of penetrating 
the interior of the brush. Consequently, for p* > 0.5, the heights of brushes formed in the systems 
with and without the presence of the fluid, become similar and when p* still increases the values 
of (/i(p*))/cr nearly flow together. The last statement is true if the segment-segment, segment-fluid 
and fluid-fluid interactions in the systems characterized by different values of and e^, are identical 
(or similar). However, when those interactions are different, then for high values of {h{p*)/a), the 
differences in the brush height in different systems remain significant, cf. figure H(b). 




Figure 6. The dependence of (h) /a on the surface density of tethered 80-mers. The nomenclature 
of the lines is as follows: T* = 1 - solid line with open circles, T* = 2 - dashed line with open 
squares, T* — 2.2 - double dash-dotted line with stars, T = 2.75 - dash-dotted line, T* = 3 - 
dashed line with triangles and T* = 4 - solid line with diamonds. The calculations are for 
= 12, and £cc ~ £■ Dotted straight line has the slope 0.35. 

The results presented in figures [4] and [5] have been computed for rather short chains composed 
of 16-mers. We now consider the case of 80-mers. Figure [6] shows the functions {h{p*)/a) evaluated 
at different temperatures for the system characterized by Ecc = £ and = 12. Similarly to the case 
of M = 16 at low surface densities of grafted chains, we observe here that the heights of the brushes 
remain almost constant. For high values of p*, the curves at different temperatures almost coincide 
and all of them can be approximated by a straight line with the slope of 0.35, in accordance with 
the theory by Alexander and de Gennes. However, we observe here an unexpected effect. Namely, 
with an increase of the temperature, the average height of the low surface density brush decreases. 
This decrease is small when the temperature rises from T* = 1 up to T = 2.2, but then the 
changes become larger, especially within the temperature range of [2.75, 3]. We carefully inspected 
that region of temperatures and found that the occurring changes are fast, but continuous, i.e., the 
plot of the thermodynamic potential Y vs T* at a constant p* does not exhibit any discontinuities. 

In order to explain the observed temperature behavior we have inspected the density profiles. 
In figure [7] we display the segment density profiles of the segment j = 2, 40 and 80 and the total 
segment density profiles. The calculations were carried out for the system presented in figure [6] 
at p* = 0.0001 and at T* = 2.2 (soUd Unes) and at T* = 3 (dashed lines). Additionally we have 
displayed here the results evaluated at the same surface density of the chains and at T* — 3, but 
assuming that the fiuid-solid energy is higher, = 18 (dotted lines). The temperature T* = 3 
is well below the bulk critical temperatures of the fluid comprising 80-mers, cf. [SS] and thus the 
chains exhibit a tendency to condensate, even within a single chain. At lower temperature T* = 2.2 
the adsorption of the fluid at the wall is stronger than at T* =3. The accumulation of the fluid 
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Figure 7. Segment density profiles psj{z) for j — 2 (a), j — 40 (b), j — 80 (c) and the total 
segment density profiles (d). The calculations were carried out for the systems from figure [o] 
at p* = 0.0001 and at T* = 2.2 (solid lines) and at T* = 3 (dashed lines). Dotted lines show 
the results evaluated at the same surface density of the chains and at T"* = 3, assuming that x 
el = 18. 



particles at the wall causes the situation wherein not all the segments of the chain "condensate", 
i.e., the initial segments do not join the "cloud" formed by the inner and outer segments and the 
tethered chains assume a characteristic mushroom-like shape. Indeed, the profiles for j = 40 and for 
j = 80 displayed in figure[7]lie very close to each other and the total segment density profile exhibits 
a depletion for 0.5 < z < 2.5 (this depletion corresponds to the mushroom's "pedicle"). When the 
temperature is raised up to T* = 3, a part of the fiuid molecules is removed from the surface 
region, making more room for the segments of the chains. Consequently, the chains are capable 
of assuming the configuration almost parallel to the wall and then they "condensate" into a thin 
(in comparison with the mushroom configuration) film covering the wall (pancake configuration) . 
However, when the wall attracts the fiuid particles more strongly, e^^ ~ 18, the formation of a thin 
film of chains at the wall is inhibited by the presence of the fiuid particles and the brush assumes 
a mushroom structure, even at T* =3. The phenomenon of the change of the brush structure is 
connected with an interplay between adsorption of the fiuid (the formation of layered structure by 
the fiuid molecules) and the tendency of chains to condensate. In the case of systems under study, 
this change is a continuous process, but we cannot exclude a situation when this transformation 
will occur as a first-order transition. 

The relations describing the changes in the brush height with the changes of the number of 
segments and with the grafting density, known as scaling relations llSffTfl |3H |35J |52] were de- 
veloped using theoretical description based on coarse-grained models. Therefore, it is not surprising 
that these relations agree well with the coarse-grained simulations (e.g., Monte Carlo simulations 
that are based on the self-avoiding random walk model of chains). However, when microscopic mod- 
elling is applied in computer simulations, several deviations between computer simulation data and 
scaling theory arise, see, for example |26) . 

In our calculations we observe that when the strength of the segment-segment interactions 
decreases, the scaling relations are approximately satisfied within the region of stretched chains. 
However, more importantly, we have found two unexpected effects. The first one is a possibility of 
the existence of a minimum on the curve describing the brush height versus the grafting density, 
while the second one is a decrease of the brush height with an increase of temperature (or an increase 
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of the brush height with a decrease of temperature) . Both effects have been observed at rather 
low grafting densities. The second effect results, in our opinion, from temperature dependence of 
adsorption of solvent molecules. In general, the effect of the adsorption of solvent on the surface on 
the brush height has been omitted in previous studies. To our best knowledge, there are no relevant 
computer simulation data that would confirm the existence of the two unusual effects mentioned 
above. Such simulations require long-lasting calculations and the question whether our observations 
are artifacts of the theory or they do really exist is currently under study in our laboratory. 
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SmIhh y CTpyKTypi 3B'5i3aHi/ix /iaHi4K)roBi/ix MO/ieKy/i ariflHO 
nepeflGaneHb Merofly 4jyHKi4ioHa/iy rycTi/iHi/i 

M. BopoBKcP, A. naxpnKee^, 0. nisid^l, C. CoKO/iOBCK|il 

Biflfliji MOfle/iKDBaHHa cfjisMKO-xiMinHi/ix npoi4eciB, yHisepcuTeT MapiT Kropi-CK/ioflOBCbKoT, 
20031 JlroGjiiH, Dojibina 

iHCTHxyT xiMifyrAIVI, 04360 MexiKO, IVleKCHKa 

Mm BMKopucTOByeMO oflHH is BapiaHxiB xeopiT 4>yHKL4ioHa/iy rycTi/iHH fl/ia BUBHeHHa 3MiH y Bucori 3B'a3aHoro 

LUapy /iaHL4tOriB, yXBOpeHHX s'eflHaHUMI/l c4)epMHHMMM CerMeHjaMU, is BMiHHUMH flOB>KMHOKD Ta noBepxHeBOKD 

rycTHHOto iiaHL4KiriB. fljia MOflejii, y jiKiCi BsaeMOfliT Mi>K cerMeHxaMU xa MO/ieKy/iaMi/i poBHUHHWKa e xaKi/iMi/i 
caMi/iMH 51K m\>k MOJieKyjiaMW po3HMHHi/iKa, HaMM BM5iB/ieHO flBa ecfjeKxi/i, aKi He cnocxepira/ii/icn y nonepeflHix 
flocjiifl>KeHH!RX. 3a neBHMX yMOB i npi/i hmbbkhx noBepxHeBi/ix KOHLieHxpaL4iiRX jiaHL4ioriB Bi/icoxa 3B'n3aHO- 
ro Luapy MO>Ke flocnraxM MiHiMyMy. KpiM xoro, fl/i5i flenKnx cucxew cnocxepiraexbciq mo npw 35ijibLueHHi 
xeMnepaxypw Bucoxa Luapy jiaHL4K)riB MO>Ke 3MeHiuyBaxMC5i. 

K/iiOHOBi c/iosa: miTKa, afl,cop6i4ia, reopifi (pyHKLiioHa/iy rycTHHH, MacwTaGyBaHHa 
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